Quantum quenches, dynamical transitions and off-equilibrium quantum criticality 
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Several mean-field computations have revealed the existence of an out of equilibrium dynamical transition 
induced by quantum quenching an isolated system starting from its symmetry broken phase. In this work 
we focus on the quantum (f> 4 A^-component field theory. By taking into account dynamical fluctuations at the 
Hartree-Fock level, corresponding to the leading order of the 1 /N expansion, we derive the critical properties of 
the dynamical transition beyond mean-field theory (including at finite temperature). We find diverging time and 
length-scales, dynamic scaling and aging. Finally, we unveil a relationship with coarsening, an off-equilibrium 
dynamical regime that can be induced by quenching from the symmetric toward the symmetry broken phase. 



Out of equilibrium quantum dynamics of isolated systems 
is a fundamental research topic which has recently become ac- 
cessible to experimental investigations by trapping ultra-cold 
atoms in optical lattices QjJ] . After the pioneering work |2|], in 
which the Mott insulator-superfluid quantum phase transition 
was observed, the field has boomed with a lot of studies, in 
particular on the so called quantum quenches. These proto- 
cols, consisting in a sudden change of an interaction parame- 
ter (for example using Feschbach resonances), bring a system 
initially in the ground state far from equilibrium. 
Out of equilibrium quantum dynamics is a very broad field. 
One of the main fascinating questions is whether, and to what 
extent, there exist universal phenomena generalizing the ones 
found for equilibrium systems. The quantum Kibble-Zurek 
mechanism, describing the production of defects occuring 
during ramps accross a quantum critical point |3J], is an ex- 
ample of such universal properties. The main topic of this ar- 
ticle is another candidate for universal behavior originally dis- 
covered in the Hubbard model [4] and later found in a large 
variety of quantum systems at the mean-field level j5)-0]. It 
consists in a dynamical transition out of equilibrium occur- 
ring after a quantum quench. Its main features are that long 
time averages display a singular behavior and the order pa- 
rameter vanishes when the final coupling Uf, reached after 
the quench, approaches a critical value Ui. 
Attempts to go beyond mean-field theory in the Hubbard 
model showed that fluctuations play an increasingly impor- 
tant role approaching the transition |^,@]. A full analysis, 
however, is still lacking. Moreover, even though it is recog- 
nized that some physical observables are singular at Ui, the 
critical nature of the transition remains to be found yet. Ac- 
tually, it is not known whether there is a diverging correlation 
length-scale at the transition nor whether some kind of crit- 
ical dynamics scaling takes place. In this work we provide 
answers to these open questions by going beyond mean-field 
theory and taking into account some dynamical fluctuations. 
In order to do that, we shall focus on the </> 4 A^-components 
quantum field theory and retain in the self-consistent 1 /N ex- 
pansion the leading contributions in the large N limit. An un- 
expected and interesting result of our analysis is that the crit- 
ical out of equilibrium dynamics occurring at the dynamical 
transition coincides with the one induced by quenches from 
the unbroken symmetry phase toward the broken symmetry 



one, a situation similar to the one leading to coarsening dy- 
namics in classical system 11011 . 

The model we focus on consists in an N component real scalar 
field interacting via a quartic term in three dimensions. It 
was studied thoroughly at equilibrium, since, depending on 
the value of N, it belongs to the same universality class as 
many physical systems such as superfiuids and ferromagnets 
flUl . The corresponding Lagrangian reads 1 121: 
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At equilibrium, this model has a quantum phase transition be- 
tween a phase with spontaneous symmetry breaking in which 
{$) is aligned along a certain direction for ro < rZ and a para- 
magnetic phase, 0) — 0, for r > r^. The critical "mass" 
ri is negative, due to the enhancement of the effective mass 
because of fluctuations. It was shown in [7] that this model 
displays, at the mean-field level, a dynamical transition due 
to quantum quenches in the mass ro (other regimes were pre- 
viously studied in 11311 ). In the following, with the aim of 
analyzing the effect of fluctuations on the dynamical transi- 
tion we retain in the two-particle irreducible/Baym-Kadanoff 
expansion of the self-energy the leading order contribution in 
I IN, which corresponds to the the dynamical Hartree-Fock 
approximation [14]. The initial condition for the dynamics is 
the ordered ground state before the quench (finite temperature 
initial conditions will be considered later). Without loss of 
generality we focus on the case where the average field, (p t , is 
aligned along the first component: <p n t = 6 nt i<f>, = (0*,}. Note 
that by symmetry the average field remains uniform for t > 
and only the diagonal terms n' — n of the connected Keldysh 
correlation functions, G™', = (f)"' f }) - <f>"<p"' , are nonzero. 
The time-dependent Dyson equations governing the evolution 
of the system after the quantum quench from ri to read: 
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FIG. 1. a) Cartoon of the dynamical transition at the mean-field level. 
From top to bottom: quench above (a), at (t) and below (b) the dy- 
namical transition, b) |r,| for a quench within the unbroken symmetry 
phase (thick blue line) and at the dynamical transition (thin red line). 
In the second case, r, decays faster than 1 jt. 



where the parallel index has been used for the n — 1 Keldysh 
correlation function and the perpendicular one for all the oth- 
ers (which are equal by symmetry). The initial condition at 
t — is given by the value of the field <p and the equal time 
(t - t' - 0) Keldysh correlation function in the ground state 
corresponding to the value of the mass ri. See the EPAPS 
for more details. Since this problem is not exactly solvable, 
we integrated numerically the equations for a large value of 
N = 10 6 (note that the average field scales as yfN). Al- 
though the dynamics of the field (p t and correlations G"" tl , look 
superficially similar to a free field evolution, the time depen- 
dence of the effective mass r, has dramatic effects as we shall 
show. 

Let us first recall the main result of mean-field theory, which 
corresponds to neglecting all the feedback of correlations on 
the dynamics of <p, in (Q} 1 1611 . The motion of the field is qual- 
itatively represented in Fig. QJ, where various quenches with 
different initial mass ri and with same final mass rl are de- 
picted. (This means that the potential V(<f>) after the quench is 
the same. The initial condition instead depends on the value of 
<p in the ground state before the quench, i.e. on ri). Above the 
transition [case (a)] the field oscillates symmetrically around 
zero and, consequently, is characterized by a zero time aver- 

— rT 

age <p = limr-»oo(l/r) J () dt<p t . Below the transition [case (b)] 
the field oscillates around one minimum of the potential and, 
hence, is characterized by a non-zero <p. In between, at the dy- 
namical transition when = r^ d ^ [case (t)] the field relaxes 
exponentially to zero, i.e. to the maximum of the potential at 
<p — 0. The phenomenology of this mean-field transition is 
identical to the one found in other mean-field models 
For example, the time averaged value of the field has a loga- 
rithmic singularity at the dynamical transition: <p oc 1 / In |A|, 
where A is the relative distance to the dynamical critical point: 



(5) 



Our goal is to determine the impact of fluctuations at first or- 
der in 1 IN on this scenario. The numerical analysis of the 
evolution eqs.(Q]|4]i shows that the system always reaches a 
steady state at long times 01711 . This is the first difference 
with respect to mean-field theory, in which oscillations instead 
persist even at long times. We show in Fig. QJ, as an exam- 



FIG. 2. a) Long time averages, <pl ^JTi and T,, as a function of the 
relative distance to the critical point A (in %). b) Critical length 
versus A the distance from the dynamical transition. Notice that 
despite the different definitions below and above the transition, f 
diverges as d/ VA on both sides of the transition, with d a and db two 
different constants. 



pie, the evolution of the mass for two different quenches: we 
find that oscillations are damped and r, converges toward an 
asymptotic value. Similar results are found for the field. By 
studying quenches for several values of the final and initial 
mass, we find that the dynamical transition continues to take 
place, as it was already mentioned in contexts related to cos- 
mology 111 811 . In the following, we study its critical features. 
Like in mean-field theory, the transition happens for quenches 
within the regime of broken symmetry: ri < r l — > ri < r[j 
and corresponds to a singularity in the asymptotic value (or 
equivalently the time averaged value) of the field. We show 
in Fig. |2k and the average mass as a function of A. Be- 
low the transition, the field relaxes to a nonzero asymptotic 
value and r t vanishes. Above the transition, the field relaxes 
to zero, whereas the mass converges to a positive value. The 
critical behavior is different from the mean-field one, since in- 
stead of a logarithmic singularity the average field vanishes as 
<p ~ |A| 1/4 approaching the transition from below (A — > + ), 
whereas the asymptotic value of r, vanishes as A for A — > 0~ 
fl^l . After having established the existence of a critical point 
let us now study its properties, i.e. focus on the physical be- 
havior after quenches right at A = 0. We find that the dynam- 
ics is divided in two stages. First, the field relaxes to zero on 
a timescale T smaller than the one characterizing the evolu- 
tion of \r,\. In the second stage, Gp tt increases exponentially, 

as G p ooe 2 ^~ p2 ~ r '', for all momenta below a cutoff A 2 = \r t= o\. 
This leads to a growth of the effective mass r,, which even- 
tually stabilizes around zero, with a slow, oscillating, power 
law decay shown in Fig.QJ). This in turn stabilizes the growth 
of Gp„. At large times, the low momentum modes enter a 
remarkable two-times dynamic scaling regime: 



cos pf 1 + 



(6) 
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with a dynamical exponent z = 1 and A a nonuniversal con- 
stant. The parallel mode G^ follows the same scaling law. 
The real space counterpart of eq. © reads Gf tt , ~ ^&(\r\ - 
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FIG. 3. a) Equal-time correlations \G plt \ as a function of p for 
t = (1000,2000,4000) in a log-log scale. Notice the divergence of 
correlations below a cutoff scale p < A ^ 0.2. b) Rescaled equal- 
time correlations \G pn /t 2 \ as a function of pt for the same data, y-axis 
in log scale. All data collapse on the scaling law {7} drawn in black. 




FIG. 4. a) Qualitative interpretation of the correlations in real-space 
in terms of a common virtual emitter in the past. G,„< vanishes in the 
dashed areas, where it is out of causal reach of virtual emitters, b) 
Rescaled two-times correlation function r G„ t i as a function of r/t 
for t/f = 1.2. All data collapse on a step function as t increases, with 
finite size effects on scale A -1 . 



(t - t'))&(t + 1' — \r\). The existence of the scaling variable f'/f 
means that the system remains always out of equilibrium: it 
is not characterized by any intrinsic time-scale besides its age 
after the quench, a phenomenon called aging [20]. 
The scaling © and (0 is demonstrated in Fig. [3] for equal- 
time correlations (in Fourier space) and for t + t' (in real 
space) in Fig.|4j). An explanation for the form of the scaling 
function can be found analyzing quenches in a free field the- 
ory where the final mass is r^ — 0. Indeed, by generalizing the 
result of ll2lll for a sudden quenches in a. free field theory we 
find the following expression for the real space two-times cor- 
relations in the continuum limit (using the notation o^ p - p 2 ): 

. , C d 3 p e'P " , , \ 

G rtf = r o — T ( C0S (<<W - f )) - cos(w p (f + t ))) . 

This is just the Fourier transform of eq. (O and 0. It's im- 
portant to realize that, contrary to the free field theory case, 
now the vanishing mass is dynamically generated by interac- 
tions. The functional form of the decrease of the mass at long 
times can be obtained, plugging the dynamical scaling form of 
the propagator into Calling A a high momentum physical 
cutoff, we find: 



r ~/ + A f dp C- 

f r A a 

= r^ + J^dp - cos(2/>0) = - 
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where to establish the last identity we have used that the con- 
stant contributions cancel since at the transition the theory is 
asymptotically massless. By taking into account sub-leading 
corrections to the dynamic scaling form of the propagator one 
can show that the mass decays even faster than 1 It 1118. 221. 
as indeed we find numerically, see Fig.QJ). Note that the map- 
ping to a free field theory, valid at large times, is also useful to 
interpret the form of the two-times scaling found previously. 
One has to use that excitations propagate at fixed speed 112 111 
and that in the limit of a large number of excitations the fields 
become classical. Then, according to the Huygens-Fresnel 
principle, plane wave propagation can be interpreted, in three 



dimensions, in terms of a continuum of virtual emitters. This 
is illustrated in Fig. [4^: between the origin and a point at a 
distance r, correlations G rt t' at successive times ?' and t are 
nonzero only provided there is a virtual emitter in the past, 
susceptible to reach the two points at times t' and t respec- 
tively. Notice that this effect includes the usual light-cone ef- 
fect found in various systems 12311 . but that the two-time 
scaling is really a new feature, due to the critical nature of all 
effective excitations. Away from the dynamical transition, we 
still observe the light cone effect but dynamic scaling does not 
hold any longer. 

We now analyze how the critical behavior emerge approach- 
ing the transition. Note that in this case there are two regimes: 
First, an out of equilibrium transient that persists for a time- 
scale t*,. In this regime, corresponding to times t such that 
r*, » t » A , the dynamical scaling © remains valid 
(on both sides of the transition) and, hence, the characteristic 
time-scale is the age of the system itself and the characteristic 
scale for the momentum is the inverse of that. In the second 
regime, corresponding to t ~ r* d , the system reaches a steady 
state in which the Keldysh correlation function becomes time- 
translation invariant. The relaxation time-scale to the steady 
state, r*,, diverges approaching the dynamical transition. Nu- 
merically we found r*, ~ 1/|A|'/ 2 . 

In the stationary regime the transverse correlation function be- 
comes time-translation invariant and has a scaling form: 



t-f 



(9) 



the low momentum behavior is critical, e.g. Gp tt ~ 1/p 2 , un- 
til values of p of the order of 1 /£* are reached, accordingly 
F(x,y) — > x 2 f(y) for x — > 0. More details on the scaling func- 
tion can be found in lfl9ll . Both r* and diverge as 1/|A|'^ 2 
approaching the transition 12411 . The fact that they are char- 
acterized by the same critical exponent is in agreement with 
the unit value of the dynamical exponent z found previously. 
The similar divergence of the de-correlation time in the steady 
state, t*, and the relaxation time toward the steady state, t*,, 
can be understood assuming that that there is no intermedi- 
ate regime. Indeed, if the out of equilibrium evolution stops 
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FIG. 5. Quench phase diagram: long time typical dynamics after a 
quench rl — > for A = 1. DT: Dynamical transition, OESB: Off- 
equilibrium symmetry breaking, R: Relaxation on large times to a 
non critical state. Error bars are smaller than item size. The exact 
position of transition lines depends on non-universal features, such 
as the interaction strength A and the cutoff A. 

when the typical momentum scale during aging, which is pro- 
portional to the time t elapsed after the quench, reaches the 
steady state value l/£*, then one finds r*, ~ l/£* ~ r*. Note 
that the asymptotic value of the effective mass r, is not directly 
related to f*. The latter is determined by studying the low mo- 
mentum properties of Gj;,, whereas the former is relevant only 
for the dependence in t — f. The diverging length is shown 
in Fig. [2b, its divergence is a power law f (A) ~ 1/A 1 ^ 2 , as 
shown in the supplementary material 

A natural question is to what extent starting from the ground 
state is important to induce the dynamical transition, espe- 
cially in regards to experiments where reaching very low tem- 
perature is one of the most challenging problem. We have ad- 
dressed this issue by considering quantum quenches in which 
for t < the system is at equilibrium at temperature T (but 
still isolated) and the value of the mass is changed suddenly 
at t — + 111 911 . We have found exactly the same critical be- 
havior provided the initial state is still in the broken symmetry 
phase. Non-universal features, such as the position of the dy- 
namical transition rl instead are different. By increasing the 
temperature at fixed r' Q one finds that the value of the critical 
mass approaches rL; they become equal when T reaches the 
value corresponding to the thermal equilibrium phase transi- 
tion. For higher temperatures the dynamical transition does 
not exists any longer. 

Let us now turn to an apparently unrelated problem: quan- 
tum quenches starting from a symmetric ground state, rl > r ( c } , 
toward values of the mass at which the system would be or- 
dered at equilibrium, ri < ri. This problem has been stud- 
ied in cosmology and in statistical physics; it is referred to 
as spinodal decomposition ifioi l25ll . Physically, one expects 
that the system globally remains in a symmetric state but lo- 
cally, on length-scales and time-scales that increase with time, 
it breaks the symmetry. Since the average field remains zero 
for all times 12611 and <p is the only dynamical quantity ana- 
lyzed at the mean-field level, the latter method is useless to 
study these quantum quenches. The growth of local order is 
visible at the level of correlations, which requires to go be- 
yond mean-field theory. In the following we briefly present 



our results obtained at the leading order in 1 /N. 
The initial conditions for quenches from the unbroken to the 
broken symmetry phase correspond to <p, = (the initial state 
is symmetric) and negative masses. These are qualitatively 
similar to those of a quench at the dynamical transition after 
the time T defined above. Indeed, it turns out that the subse- 
quent out of equilibrium dynamics is the same. In particular, 
the effective mass vanishes asymptotically, and the two-time 
correlations scale like (|7). Thus, at least at the leading or- 
der in 1 IN, the dynamical transition is characterized by the 
same critical properties of coarsening dynamics. A complete 
quench phase diagram is shown in Fig. [5^, summarizing all 
possible quenches r' Q — > r^ . When the initial field is nonzero, 
ri. < r[j, the system relaxes to a steady state on both sides of 
the dynamical transition, either to a state of positive field (f> or 
of positive mass T,. The correlations follow the scaling form 
© on the dynamical transition (DT) and in the whole region 
(OESB) of quenches from the symmetric phase to the broken 
symmetry phase. 

In conclusion, by going beyond mean-field theory and tak- 
ing into account fluctuations at the leading order in 1 jN we 
have shown that the dynamical transition induced by quantum 
quenches displays bona fide critical properties, in particular 
diverging time and lenth-scales, dynamical scaling and aging. 
Moreover we have found that this aging dynamics is identi- 
cal to the coarsening dynamics following quantum quenches 
from the symmetry broken to the symmetry unbroken phase. 
Clearly, it would be interesting to extend our analysis to mod- 
els directly relevant for experiments, such as the Bose Hub- 
bard model, and to the next leading order contribution in 
l/N, where several different equal-time scaling fixed points 
were found and long-time thermalization, at least for small 
quenches, is guaranteed [27]. This would provide a crucial 
test to establish the nature of the dynamical transition. We 
envision two different possibilities: either the transition per- 
sists when dynamical fluctuations inducing quantum noise are 
included; in this case it is likely related to the thermal phase 
transition crossed by increasing the strength of the quench, 
or, more interestingly, it becomes a cross-over related to a dy- 
namical phenomenon that takes place while the system is tran- 
siently trapped in a off-equilibrium region of the configuration 
space. 
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